Graphene integer quantum Hall effect in the ferromagnetic and paramagnetic regimes 
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Starting from the graphene lattice tight-binding Hamiltonian with an on-site U and long-range Coulomb 
repulsion, we derive an interacting continuum Dirac theory governing the low-energy behavior of graphene in 
an applied magnetic field. Initially, we consider a clean graphene system within this effective theory and explore 
integer quantum Hall ferromagnetism stabilized by exchange from the long-range Coulomb repulsion. We study 
in detail the ground state and excitations at ;/ = and u = ±1, taking into account small symmetry-breaking 
terms that arise from the lattice-scale interactions, and also explore the ground states selected a.\.v = ±3, ±4, 
and ±5. We argue that the ferromagnetic regime may not yet be realized in current experimental samples, 
which at the above filling factors perhaps remain paramagnetic due to strong disorder. In an attempt to access 
the latter regime where the role of exchange is strongly suppressed by disorder, we apply Hartree theory to study 
the effects of interactions. Here, we find that Zeeman splitting together with symmetry-breaking interactions 
can in principle produce integer quantum Hall states in a paramagnetic system at z/ = 0, ±1 and ±4, but not 
at I' = ±3 or ±5, consistent with recent experiments in high magnetic fields. We make predictions for the 
activation energies in these quantum Hall states which will be useful for determining their true origin. 

PACS numbers: 73.43.-f, 71.10.-w, 71.10.Pm 
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I. INTRODUCTION 

Recent experimental advances have made possible the iso- 
lation of high quality two-dimensional graphene sheetsjl thus 
opening up a new arena for exploring quantum Hall physics. 
As illustrated in Fig. ^ graphene is a collection of carbon 
atoms arranged on a honeycomb lattice. At half-filling and 
in the absence of a magnetic field, graphene is semi-metallic 
with a Fermi surface consisting of two distinct nodes resid- 
ing at the Brillouin zone corners. Low-energy excitations 
about the two nodes are well characterized by two "flavors" 
of Dirac fermions whose energies obey a linear dispersion re- 
lation, i?(k) = ±fi,i;|k|, where for graphene v « lO^m/s.^-^ 
Applying a magnetic field rearranges the spectrum into Lan- 
dau levels, each of which is approximately four-fold degen- 
erate near the Fermi level owing to the presence of flavor 
and spin. The underlying Dirac structure gives rise to an un- 
conventional quantum Hall effect which has been observed 
experimentallyii^ii^ and explored theoretically from numer- 
ous perspectivesi^i'^i^'^i^°i"i'^i^^i^^i'^i'^ 

Initial experiments on the graphene quantum Hall effect re- 
ported quantized Hall plateaus at cr^tj = ve^ jh = 4(n + 
\/2)e'^ /h, where n is an integer. ^••' The appearance of these 
integer quantum Hall states is quite natural from a non- 
interacting, single-particle perspective, given the approxi- 
mate four-fold Landau level degeneracy. There are, however, 
at least two mechanisms by which additional integer quan- 
tum Hall states in graphene may be induced. First, such 
states can in principle be stabilized by explicit symmetry- 
breaking fields and/or interactions, the simplest mechanism 
being Zeeman splitting. In this situation the system would 
be appropriately characterized as a "quantum Hall param- 
agnet". Second, provided disorder is sufficiently weak the 
poorly screened Coulomb interactions in graphene can dra- 
matically modify the single-particle picture, in analogy with 
the well studied "quantum Hall ferromagnetism" in GaAs 



heterostructures.'^ In the latter system, for instance, exchange 
interactions strongly stabilize the integer quantum Hall state at 
V = 1, even in the limit of vanishing Zeeman coupling — that 
is, the ground state is a ferromagnet. Generally speaking, this 
mechanism is actually far more common than the first since 
the Coulomb energy scale is typically much larger than the 
energy scale associated with symmetry-breaking terms such 
as the Zeeman energy. Additional exchange-driven integer 
quantum Hall states can similarly be expected to appear in 
a sufficiently clean graphene system. In fact, quantum Hall 
ferromagnetism in graphene promises to be even richer than 
in GaAs due to the additional flavor degree of freedom, which 
as we will discuss offers the interesting possibility of lattice- 
scale order coexisting with the integer quantum Hall effect. 
The problem is thus more analogous to quantum Hall bilayer 
and Si MOSFET physics,'** where layer and valley indices 
play a role similar to flavor in graphene. 

Interestingly, more recent experiments utilizing higher 
magnetic fields have in fact resolved additional quantized Hall 
plateaus at filling factors v — Q, ±1, and ±4.^ The activation 
energy measured at i/ = ±4 was found, somewhat surpris- 
ingly, to be dominated by the single-particle Zeeman energy, 
suggesting that at these filling factors the system resides in the 
paramagnetic regime due to strong disorder The origin of the 
quantum Hall states at t^ = and j/ — ±1 is unclear at present, 
and presents an interesting puzzle. In a very recent paper No- 
mura and MacDonald suggest that these states are due to the 
onset of quantum Hall ferromagnetism.'^ Given the apparent 
spin-splitting origin of the states at j/ = ±4, however, it is per- 
haps worth exploring an alternative scenario where the quan- 
tum Hall states at i^ = and v = ±1 are similarly due to 
explicit symmetry breaking, rather than Coulomb exchange. 

The purpose of the present paper is twofold. In the first 
part, we carry out a detailed exploration of quantum Hall fer- 
romagnetism, focusing on a clean graphene sample. Follow- 
ing this analysis, we consider a dirty system and ask whether 



symmetry-breaking terms can stabilize integer quantum Hall 
states in the paramagnetic regime. To this end, we start 
with the lattice Hamiltonian including on-site U and long- 
range Coulomb interactions, and derive a continuum interact- 
ing Dirac formulation suitable for studying these phenomena. 
In addition to the usual long-range Coulomb repulsion, we re- 
tain shorter-range terms that arise from the microscopic inter- 
actions and explicitly break the flavor degeneracy exhibited in 
a non-interacting theory. Such symmetry-breaking terms have 
not been taken into account previously, and are important for 
analyzing the physics both in the ferromagnetic and param- 
agnetic limits. In the ferromagnetic regime, we analyze the 
ground state and excitations within this theory at integer fill- 
ing factors in the lowest two Landau levels. At filling factors 
ly = ±1, we show that interactions favor "easy-axis" flavor 
polarization, leading to a ground state exhibiting charge den- 
sity wave order on the lattice scale. The i^ = ground state is 
found to depend on the strength of interactions relative to the 
Zeeman coupling, and will either be a uniform spin-polarized 
state or a spin-singlet exhibiting the same lattice-scale struc- 
ture as at !^ — ±1. We calculate the spin- wave and particle- 
hole excitation energies at these filling factors, as well as iden- 
tify the relevant skyrmions which are expected to set the trans- 
port gap. Interestingly, interactions favor "easy-plane" flavor 
polarization aiv — ±3 and ±5, leading to a ground state that 
spontaneously breaks U(l) flavor symmetry. Thus, a finite- 
temperature Kosterlitz-Thouless transition can be expected at 
these filUng factors. Finally, a uniform, spin-polarized ground 
state is expected aX. v = ±4. We establish that skyrmions 
continue to provide the minimum-energy charge excitations 
?A V = ±3, ±4, and ±5 as well, at least in the absence of 
anisotropy terms. 

To explore the paramagnetic regime in a dirty system, we 
use the expectation that the effects of exchange should be 
strongly suppressed by disorder and incorporate interactions 
within Hartree theory. Here, we argue that in principle in- 
teractions can give rise to quantum Hall states siv — 0, ±1, 
and ±4, but not at ly — ±3 or ±5, which is consistent with 
the recent high-field experiments. We provide estimates for 
the activation energies in these additional quantum Hall states, 
which will be useful for determining experimentally whether 
the 1/ = and ±1 states originate from explicit symmetry 
breaking or quantum Hall ferromagnetism. 



II. CONTINUUM INTERACTING THEORY 

We start with the zero-field lattice Hamiltonian for 
graphene written as 



H — Ht + Hu + ifcoui, 



(1) 



where Ht describes hopping of electrons across nearest- 
neighbor honeycomb sites, 

Ht = -tY. ^ [cLcax' +H.C.], (2) 

(xx') a=T,J. 



and Hij and -ffcoui contain the on-site repulsion and long- 
range Coulomb interaction, respectively: 
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(3) 
(4) 
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Here n^ = n^^ + n^x is the electron number operator and 
S(x) = ^Cax^a^cpx IS the usual spin operator with cr a 
vector of Pauli matrices. The Coulomb potential is V{x) — 
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Tne M' ^^^ ^ ^^ appropriately chosen dielectric constant (see 
below). We will introduce the magnetic field upon entering 
the continuum, which is justified since for experimental field 
ranges the magnetic length ts — \/ft/(ei?) is much larger 
than the lattice spacing ap. 

As mentioned above, with one electron per site, the honey- 
comb band structure exhibits two gapless Dirac points at the 
Fermi energy, occurring at wave vectors ±Q — ±(47r/3, 0). 
Focusing on the linearly dispersing excitations near the two 
nodes, the Fourier-transformed lattice fermion operators may 
be conveniently expanded in terms of two flavors of contin- 
uum Dirac fermion fields (denoted R and L) as follows. 



Caq+Qa ~ 7V'afla(q) 



(5) 
(6) 



Here and below we reserve indices a, (3 for spin, A, B for 
flavor, and a, b for the honeycomb sublattice. Moreover, 
CT^o, T^^, and 77;^^ denote Pauli matrices that contract with 
the spin, flavor, and sublattice indices, respectively. (Note 
that according to Eqs. (|5} and (|5} V'i?! and '4>l2 correspond 
to sublattice 1, while ':J}r2 and ipLi correspond to sublat- 
tice 2.) We will use the convention that suppressed in- 
dices on the fields are implicitly summed {i.e., ip'^ip = 
J2aAa '^l.Aa'^o'Aa)- The normalization on the right side of 
Eqs. (|5} and (|6} is chosen to be 7 = \/2/{3^^^L) so that the 
continuum fields satisfy the canonical anticommutation rela- 
tions {■IpaAaicO^'^'lBbi^')} = Sa/sSABSabi^Tr)^ 5^ {q ~ q'). 

Real-space Dirac fields are defined according to 



^QAa(q) 



d^xe 



^IpaAai'X.) 



(7) 



and satisfy {V'aAa(x), V^Bh(x')} = SaisSAsSabS^ix- x'). 
Note that consistency requires ip{x) to be a slowly varying 
field, or equivalently that V'(q) is nonzero only within some 
momentum- space cutoff A. The lattice electron number op- 
erator Tixa for sublattice a can now be expressed in terms of 
real-space Dirac fields as follows. 



nxi 



nx2 
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PRi + PL2 + e^'^^-^jl + e-2'Q- J+] (8) 



[PR2 + PL 



^ _ g2iQ.x j-^ _ g-2.Q.x j_^9) 



where pAa = Ea=Ta ^^Aa^^aAa Contain uniform pieces of 
the density and J+ = 'ip\ii4'L2 and J_ = Wr2^li contain 




FIG. 1: Schematic lattice-scale order in the u = —1 integer quan- 
tum Hall state. "1" and "2" label the two sublattices of the honey- 
comb. Electrons in the highest-occupied Landau level reside only on 
sublattice 1 (for instance), and undergo circulating currents around 
second-neighbor plaquettes oriented along the arrows. 



oscillatory components of the density. The spin operators on 
each sublattice take on similar forms, with ^cr inserted be- 
tween each V'^ and ijj above. It follows that in the continuum, 
the uniform part of the total density is ptot = ip^ip^ while the 
uniform spin density is Stot = ■^ip'^crtp. 

Using the continuum expansion above and turning on the 
external magnetic field, we arrive at the desired continuum 
Hamiltonian, which we write as 



Ti. = Ti-o + Wi. 



(10) 



Here Hq denotes part of the Hamiltonian that is invariant un- 
der global 5'C/(4) = SU{2)spin X S'[/(2)fiavor rotations. 



+ i /"d2xd2x'ptot(x)^(x-x')ptot(x'), (11) 

with Dj — dj — i{e/h)Aj. Hi encodes the remaining 
anisotropy terms which break this SU{4:) symmetry down to 

C/(l)spin X [C/(l) X Zalflavor, 



Hi = / d^xj -g/isB-Stot + ^uo[Ptot +P; 



2 

'staff 



- -(S2ji+sL + 6S;,.i.Si2 + (l--2))] 



Y^ ui(r)pstag(x + r)pstag(x) 

r 



(12) 



where Pstag = ip'^T^ri^ilj represents the staggered electron 
density between sublattices 1 and 2 of the honeycomb and 
Sag — \wAa'^'^Aa- The sum in Eq. ( fT2t is over triangular 
lattice vectors r — me+ + ne_, where m, n are integers and 
e± = (ao/2)(±l, VS)- The first two lines of Eq. (I12> contain 



the Zeeman coupling and the continuum form of the on-site 
U, with 



Mo = y/lalU/A. 



(13) 



For concreteness we take the magnetic field B perpendicu- 
lar to the graphene plane along the +z direction throughout. 
Apart from the long-range density-density repulsion exhib- 
ited in Tio, the lattice Coulomb interaction also gives rise to 
the shorter-range terms in the last two lines of Eq. ( I12> . re- 
flecting lattice-scale physics. The vi term represents an inter- 
sublattice repulsion that reflects the smaller Coulomb energy 
cost for electrons residing on the same sublattice versus oppo- 
site sublattices. The M2 term represents the intra-sublattice re- 
pulsion between oscillating components of the density in Eqs. 
(|8} and (|9}. The coupling constants for these terms are given 
by 



vi{v) 



Vial 



U2 



[V{v + 1/V3y) - (1 - 6,^o)V{v)] (14) 

15) 



— -^V^(r)cosQT~ '■^' 



r^^O 



3 °V4^eao 



For later use we note that 



ui = ^fi(r) « —U2- 



(16) 



Some comments are worthwhile here. The above contin- 
uum formulation is expected to remain appropriate out to a 
cutoff A ^ 7r/(4ao) or so from the Dirac nodes, which con- 
stitutes an appreciable fraction of the Brillouin zone. Elim- 
ination of modes outside of this radius will renormalize the 
parameters in the theory somewhat, though for simplicity we 
have retained the bare values for the coupling constants. Fur- 
thermore, other symmetry-allowed interactions will in prin- 
ciple be generated as well; however, we shall assume these 
are subdominant. We also remark that in Eq. J12l i we have 
only retained terms arising from the lattice Coulomb repul- 
sion that are nonvanishing upon assuming a local form in the 
continuum (namely, vi and U2); other terms are expected to 
essentially average out and will thus be unimportant. Given 
that ?/; is a slowly varying field, it is tempting to approximate 
vi as a local interaction (as we have done for U2). For most 
purposes this is in fact adequate, and we shall often make this 
assumption to simplify our analysis. However, as we will dis- 
cuss below, despite appearances assuming a purely local form 
of vi leaves the SU(2) flavor symmetry unbroken in the ferro- 
magnetic ground state at filling factors v = ±1, which is why 
we have retained the finite range of this interaction. 

We now turn to screening in the Dirac theory. Since the vi 
and f 2 terms are effective only on short length scales, these 
interactions are essentially unmodified by screening, and we 
thus assume that the bare dielectric constant for graphene ap- 
pears in Eqs. (I14> and ( I15> . With air on one side of the 
graphene plane and Si02 on the other, the unscreened dielec- 
tric constant is estimated to be e « [(^+^ + ^-^)/2]^^ ~ 
1.6eo, where e+ w eo and e_ w 4eo correspond to air and 



Si02, respectively. The long-range part of the Coulomb inter- 
action in Ho is, however, expected to be (weakly) screened. 
To estimate the screening, we first note that the presence 
of the applied magnetic field is expected to become im- 
portant on length scales longer than roughly the magnetic 
length, which is much larger than the inverse cutoff A^^. 
Thus the screened dielectric constant for the zero-field case 
should not be dramatically modified by the presence of the 
field. We therefore use the zero-field screened dielectric con- 
stant for the long-range Coulomb repulsion in Eq. Jlll i. ob- 
tained within the random phase approximation'^ which yields 



fRPA 



+ e^/{8hv) «5eo- 



IV. INTEGER QUANTUM HALL FERROMAGNETISM 

In this section we focus on a clean graphene system at 
integer filling factors, and explore "quantum Hall ferromag- 
netism" using the interacting theory derived above. We will 
first analyze the n ~ Q Landau level, and then turn our at- 
tention to the n = 1 Landau level. In the last part of this 
section we briefly discuss current experimental relevance for 
our results. 



A. n = Landau level 



III. OVERVIEW OF DIRAC LANDAU LEVELS 

Before analyzing the full interacting theory, it will be use- 
ful to briefly recall the Landau level spectrum for the non- 
interacting case.*'^" Throughout the paper we implicitly work 
in the symmetric gauge. As a consequence of the Dirac struc- 
ture, the Landau level energies in the non-interacting theory 
are given by 



E 



Ui:' 



== T^ff/^sS + sign(n) V2e&^2^ 



(17) 



where n is an integer Zeeman coupling weakly spin splits 
the Landau levels, leaving a two-fold flavor degeneracy in the 
absence of interactions. The two-component wave function 
for spin a and flavor A in the ?i*'' Landau level can be written 






X" 



For n 7^ 0, the elements </)" and x" are simply wave functions 
for Schrodinger equation Landau levels with indices |n| — 1 
and \n\, respectively. In the symmetric gauge, the elements are 
related by x" = \n\^^/^sign{n)a^(l)^\ where a^ is the usual 
raising operator for Schrodinger equation Landau levels. One 
can show that the total probability weight is divided equally 
between the upper and lower elements of the wave function; in 
other words, an electron with wave function $^^ is equally 
likely to be found on either honeycomb sublattice. In contrast, 
the n = Landau level wave functions have 0"=" = 0, while 
-^n=o jg ^ Schrodinger equation lowest Landau level wave 
function. By examining Eqs. (|5) and (|5Jl we see that ^^^]f has 
weight only on sublattice 2 of the honeycomb, while $^2*^ 
has weight only on sublattice 1. This qualitative distinction 
between the n — and n ^ Landau level wave functions 
has important consequences when we incorporate interactions 
below. 

We note that this feature of the n — Landau level 
wave functions is robust against introducing spin-orbit cou- 
pling, which admits a term TCso — '^so J (P:>i.ip^ri^(J^ip to 
the Hamiltonian^'. In contrast, in the presence of this term 
the probability density for n ^ Landau level wave func- 
tions will no longer be evenly distributed between both sub- 
lattices. This correction is quite small, as the induced den- 
sity difference between the sublattices will be proportional to 
^so^bH^v) <C 1, and will hereafter be neglected. 



We will employ two simplifying assumptions in our anal- 
ysis of quantum Hall ferromagnetism in the n — Q Landau 
level. First, we make the standard approximation of ignor- 
ing Landau level mixing and project out states away from 
the 71 = Landau level. This is reasonable given that the 
spacing between the n = and n = ±1 Landau levels is 
roughly 420^/b|tJk, where B[T] is the magnetic field in 
Teslas, while the Coulomb energy scale 



£c 



47reRPA^B 



130VS[T]K 



(19) 



is around three times smaller Second, since the non- 
interacting wave functions live on only one sublattice or the 
other, we will set ipaR/Li ^ 0. With the electron kinetic 
energy quenched, TLq then becomes 



(18) Hr° ^\j d'xd2xVt„t(x)F(x - x'; 



Ptot (,X 



(20) 



The anisotropy terms in Tii are dramatically simplified upon 
projection. 



7^r° = /"d^xj-g^sB-St 



2 



'A2 



'A2J 



A=R/L 
E I'lWPstaglx + r)pstag(x)}. 



(21) 



In the projected subspace, ptot = PL2 + PR2, Stot = Sl2 + 
S7J2, and Pstag — Pl2 ~ Pr2- Note that the U2 interaction in 
Eq. ( I12> has dropped out altogether 

In the SU(4)-invariant limit with 7i"=° = 0, despite the 
four-fold degeneracy of the n = Landau level the system 
is an exchange ferromagnet at filling factors v = —1, 0, and 
+1 (corresponding to a quarter-filled, half-filled, and three- 
quarter-filled level), and exhibits a quantized Hall conduc- 
tance a^y = ve^ jh. The SU(4) spin/flavor symmetry is bro- 
ken spontaneously in the absence of anisotropy terms, though 
in accordance with the Mermin- Wagner theorem.^^ true long- 
range order can exist only at zero temperature. This spon- 
taneous symmetry breaking gives rise to three gapless Gold- 
stone modes at t/ = ±1 and four at j/ = 0. These SU(4) "spin- 
waves" have the following dispersion at small momentum^ 



£;o(q) = «(4g)', 
where the stiffness is given by 



"^iVf^^- 



(22) 



(23) 



The lowest-energy charge excitations in these integer quantum 
Hall states are known to be SU(4) skyrmionsiiSlii^ whose 
energies are determined by the stiffness k. These topological 
excitations were explored in Ref. 23. Due to the Coulomb 
repulsion the minimum-energy skyrmions are infinitely large, 
and the energy cost for creating a skyrmion/antiskyrmion pair 



E. 



sk/ask 



2k, 



(24) 



which is half that for an electron-hole pair 

The goal of the remainder of this subsection will be to 
address how this picture is modified by anisotropy terms in 
7i"^°. In particular, we will discuss the nature of the ground 
states selected by these symmetry-breaking terms and exam- 
ine their effects on excitations out of the resulting ordered 
states. 



1. Filling factors v — iil 

Consider now the v = —\ state, corresponding to a quarter- 
filled n = Landau level {v = +1 is related by particle-hole 
symmetry and will not be discussed separately). Zeeman cou- 
pling favors a ground state occupied solely by spin up elec- 
trons. Moreover, the vi term provides an "easy-axis" flavor 
symmetry that favors having (pstag) 7^ by occupying ei- 
ther all flavor R or all flavor L states. It follows that in the 
presence of interactions the v ^ —1 ground state exhibits a 
spontaneously broken Z2 flavor symmetry, and can therefore 
sustain long-range order at finite temperature. Filling, say, the 
flavor R states, the ground state can be written 



-i)-n4 



Tfl.ml 



where \vac) is the fermion vacuum and c 



(25) 



t 



adds an elec- 
tron in the 71 = Landau level with spin a, flavor A, and 
angular momentum m. Since flavors L and R correspond to 
sublattices 1 and 2 in the projected subspace, the i^ = — 1 
ground state is a spin-polarized charge density wave (CDW), 
having a larger electron density on one of the two honeycomb 
sublattices. The charge order is illustrated schematically in 
Fig. n Physically, this CDW arises because electrons can 
remain farther apart by occupying one sublattice before the 
other, thereby minimizing their Coulomb energy. We em- 
phasize, however, that only the relatively small number of 
electrons participating in the n — Landau level (roughly 
1.4 X 10^'^i?[T] electrons per hexagon) contribute to this or- 
der, which might make it difficult to observe with an STM. 



Another experimental signature of the CDW at i^ = — 1 may 
come from carbon NMR measurements. First, there is a net 
electron spin only on the higher-density sublattice, producing 
a magnetic moment per site of roughly 1.4 x 10^^B[T]ijLb 
with which the (relatively rare) carbon atoms with nuclear 
spin can interact. A second effect arises from more subtle 
lattice-scale structure that exists in the v = —\ state. Indeed, 
it follows from Eqs. (|3 and (|6j that electrons on the higher- 
density sublattice participate in currents circulating around 
second-neighbor plaquettes, oriented along the arrows in Fig. 
n This circulation in turn gives rise to a magnetic moment 
along the applied field direction that couples to atoms in the 
lower-density sublattice. Using our zero-field relation be- 
tween the lattice and continuum fields, we crudely estimate 
the magnitude of this induced moment to be 5 x 1Q^^B\T]^ib- 
These effects should induce a small splitting in the nuclear 
spin precession frequencies for carbon atoms on the two sub- 
lattices, which may be measurable. 

We turn next to excitations in the presence of anisotropy 
terms. The three Goldstone modes present in the SU(4)- 
invariant limit become gapped. These generalized "spin- 
wave" gaps contain important information regarding the sta- 
bility of the spin and CDW order in the ground state, as well 
as for the transport gap. Employing the usual lowest Landau 
level projection, the energies for these three branches can be 
obtained essentially exactly. We follow closely the approach 
in the review in Ref. 17, generalized to the SU(4) case. The 
procedure is straightforward though tedious, and here we shall 
only outline the calculation. Define 



S^"" = ^ij^'T'^a''^, 



(26) 



where /z, ly run over 0, x, y, z and t^ , a° are identity matri- 
ces. We wish to calculate the energies of S'^^{i:\)\v = —1), 
^^{q)W = -1)' and 5^(q)|i/ = -1). which turn out to 
be exact eigenstates of the projected interacting Hamiltonian. 
Here the overbar indicates a lowest Landau level projection, 
S'O^(q) creates a spin-wave with momentum q, 5^°(q) cre- 
ates a "flavor- wave", and S^^ (q) creates a mixed spin/flavor- 
wave. The energies of these excited states are conveniently 
calculated in first quantization. The Hamiltonian in momen- 
tum space becomes 



-T/n— () 

I-Lq 

-T/n— 



1 



y(k)pt„t(k)ptot(-k) 



(27) 



-gfisBSO'^Ck^O) 



1 



+ -uo I [pt„t(k)p,„t(-k)+45-"(k)5-"(-k) 



- -(S'^*(k)5^^(-k) + 50^(k)S'0*(-k))] 

- 4^t-i(r) / e*''"-5^(k)^(-k), 



(28) 



where i is summed over x, y, z. The projected first quantized 



Ptot(k) = e 



E^kO) 



operators are given explicitly by 

3 

5^(k) = ie-^(^«^-)'EPkO>^(j>''(j), (29) 
i 

where j runs over each electron in the ground state and pk(i) 
is a unitary operator satisfying 



Pk(j)Pq(j) 



e' 2 



■fcAg 



Pk+q(i) 



(30) 



bk(i),Pq(j)] = 2i,5ySin(4fcAg/2)pk+q(i) (31) 

with k f\ q = z • (k X q). Evaluating the commutators 
[Ho=° + H'i'=°, S^(q)], efc. on the ground state yields the 
corresponding excitation energies. 

At small momentum, we obtain the following energies for 
spin-wave, flavor-wave, and mixed spin/flavor-wave excita- 
tions, respectively: 



EficL) 



g^isB + naifsqf 




a-o 



uipo + Kfiisgy 



24V^B, 
Eafici) - gfiBB + Ef{q), 

where po = ^/{'^'^^'b) and the stiffnesses are 

Ka = (Uo - Ui)po + K 
Kf = UiPo + K. 



(32) 

(33) 
(34) 



(35) 
(36) 



At large momentum |q| — > cx), these excitations correspond to 
the three types of separated particle-hole excitations out of the 
ferromagnetic state (particle-hole excitations can be created 
by flipping a spin, flavor, or both). Their energies are 



£;cr(oo) = gp.BB + 2{uo -ui)po 
Ef{(X)) = 2uipo + Ak 
Eaf (oo) = gfisB + 2uipo + Ak. 



4k 



(37) 
(38) 
(39) 



The most noteworthy feature of these excitation energies is 
the small value of the gap for flavor-wave excitations, which 
is Ef{0) « 4 X 10-^{B[T]f^^K. In fact, had we ignored the 
finite-range of vi here, we would incorrectly conclude that the 
flavor-wave gap vanishes. Assuming a local form, the sub- 
lattice repulsion becomes proportional to S'^'^(x)^, which is 
simply a constant when acting on uniform density states in 
the quarter-filled Landau level. The lack of symmetry break- 
ing from such local interactions has been noted previously 
in studies of quantum Hall bilayersi^ As a result the flavor 
wave gap is down by an additional factor of ao/is from what 
one might naively anticipate. Nevertheless, long-range CDW 
order is expected to persist up to a transition temperature of 
order Kf/\n[Kf/Ef{Q)]}'' The spin order is comparatively 
more robust, as the spin-wave gap is gpsB ^ B\T]K,^ which 
for typical field ranges is around 50 times larger or so than the 
flavor wave gap. 
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The smallness of the flavor-wave gap has important impli- 
cations for skyrmions as well, which will continue to set the 
charge gap. In particular, skyrmions will be cheaper to cre- 
ate solely in flavor space, leaving the spin order intact. The 
competition between the sublattice repulsion vi, which favors 
small skyrmions, and the long-range Coulomb, which favors 
infinite skyrmions, will set a length scale for the minimum- 
energy skyrmions. However, the ratio of the flavor-wave gap 
to the Coulomb energy scale Sc is quite small, given ap- 
proximately by 3 X 10~^i3[T]. For comparison, the ratio of 
the Zeeman gap to the Coulomb energy in GaAs is roughly 
6 X \0-^ ^ B[T],'^^ which at lOT is around 60 times larger. 
Hence, corrections to the energies of flavor skyrmions arising 
from the nonzero flavor-wave gap are expected to be small. 
Approximating their energies by those for infinite skyrmions 
in the absence of a flavor- wave gap, the activation energy aris 
ing from flavor skyrmions is estimated to be Ay=_i w k/ 
These results imply that the CDW order will diminish ex- 
tremely rapidly as one tunes away from v = —1, whereas 
the spin order will be much more robust. One possible test of 
these predictions would be to measure the activation energy at 
constant perpendicular magnetic field with varying in-plane 
fields-^. The activation energy arising from flavor skyrmions, 
being independent of spin, will be insensitive to changes in 
the latter 



2. Filling factor u = 

The ground state at i^ = 0, corresponding to a half-filled 
n = Landau level, will depend on the strength of the Zee- 
man coupling and on-site U relative to the sublattice repulsion 
vi. If the former terms dominate, the ground state will be a 
spin-polarized flavor-singlet. 



o)sp=n4 



i?,m4i,™l^"^) 



(40) 



which in contrast to i^ = — 1 respects all graphene lattice sym- 
metries. If the sublattice repulsion dominates, however, the 
ground state will be a flavor-polarized spin-singlet. 



0>fp=n4fl,m4ii,ml""'=) 



(41) 



and thus also exhibit the lattice-scale CDW order illustrated 
in Fig. n Comparing the energies of these two states, 
one finds that the ground state will be spin-polarized if 
gpsB + 2uqPo > Auipo, while a flavor-polarized ground 
state emerges if gpsB + 2uqPo < 4uipo- According to 
Eqs. ([13} and ^, we have uopo ~ 0.08C/[eV]B[T]K and 
uipo w 0.4i3[T]K, where [/[eV] is the strength of the on-site 
repulsion evaluated in electron-volts. Ascertaining which sce- 
nario prevails is complicated by the rather large uncertainty in 
the value of the on-site U. Previous estimates have suggested 
U ^ 5-12eV,^'^ though here a slightly lower range of values 
may be appropriate due to the presence of the Si02 substrate. 
With g = I A, which is the smallest value reported in recent 
high-field experiments^, the critical value for U above which 



the ground state is spin-polarized is estimated to be Uc ^ 
4eV. Hence we can not conclusively determine whether the 
ground state is spin-polarized, though this regime can always 
be reached by increasing the effective Zeeman coupling via 
the introduction of an in-plane magnetic field component. Be- 
low we will discuss excitations in both regimes. 

Consider first the situation where the ground state is spin- 
polarized. Here there are four low-lying branches of excita- 
tions, generated by 5^(q) ± iS^(q) and S^(q) ± ^S^(q), 
which create spin waves by breaking flavor singlets as re- 
quired by the Pauli principle. Schematically, the former 
breaks singlets by sending {\R)\L) - \L)\R)) ^ \R)\R) (or 
\L)\L)\ while the latter sends {\R)\L) - \L)\R)) -^ \R)\L) 
(or \L)\R)). That is, S'^^(q) ± iSv^{q) creates a spin-wave 
by doubly-occupying a given flavor, while S''^^(q) ± S'^^(q) 
does not. Proceeding as in i^ = — 1, at small momentum the 
respective energies for the first and second pairs of spin waves 
are 



ETM) 

i^r±(q) 



giisB + 2uoPo - 4uipo + t^fi^Bq)'^ (42) 
giiBB + K^iieq?, (43) 



(44) 



while at |q| — > oo we have 

-B^^2±(oo) = g^^BB + 2mo/3o ~ 2-U1/90 + 4k. 

The stiffnesses K,f and n^ are the same as at ;/ = —1. Note 
that having _E^!^(0) > requires that gjisB + 2uopQ > 
4:Uipo, which is the same condition given above for realizing 
a spin-polarized ground state. 

We will assume that uq — 2ui < 0, implying that El^ (0) < 
-E2±(0). This occurs for U < lOeV, which we believe is rea- 
sonable. Due to the smaller gap for exciting spin waves by 
doubly-occupying a given flavor, it follows that the charge gap 
will be set by mixed spin/flavor-textured skyrmions. Unfortu- 
nately, the large uncertainty in the strength of the on-site U 
does not permit us to make a quantitative estimate for the acti- 
vation energy set by these skyrmions. We can at least say qual- 
itatively that as one tunes away from v = Q they should induce 
a depletion in the total spin of the system along with a simul- 
taneous revival of the sublattice CDW order How rapidly this 
occurs will of course depend on the optimal skyrmion size, 
which we are unable to determine. 

Consider now the case where the sublattice repulsion dom- 
inates, leading to a flavor-polarized ground state. There are 
again four low-lying branches of excitations, generated by 
^S^(q) ± iS^(q) and ^(q) ± ^S^(q), which create flavor 
waves by breaking spin singlets. The first pair of operators 
break singlets by sending (| T)U) - I i>l T» ^ I T)l T> 
(or I i)| i)), while the latter operators send (| T)l i) - I i 
>l T>) ^ I T>U> (or U)| t»- The energies for the first and 
second pairs of flavor waves are 

£;lP±(q) « Tg/isS + 2(2ui-uo)po + K/(^i3<z)'(45) 
£;*(q) « 2{2u^-uo)po + nf{lBqf (46) 

at small momentum, while 

-B^j-(oo) ^ =F5MsS + 2(3ui - Mo)po 

-E2±(oo) = 2(3mi — wo)po + 4k. 



Here the gap is smallest for flavor waves exciting by dou- 
bly occupying spin up states as favored by Zeeman coupling, 
so the charge gap will be set by mixed spin/flavor-textured 
skyrmions in this case as well. Away from v = {) such 
skyrmions will produce a rapid depletion of the CDW order 
and revival of spin polarization. 

Distinguishing experimentally between these two ground 
states can be achieved by measuring the activation energy 
in the presence of an in-plane magnetic field that enhances 
the effective g-factor If the ground state is spin-polarized, 
then the activation energy should increase with the in-plane 
field since the larger g-factor favors smaller (and thus higher- 
energy) skyrmions. On the other hand, if the ground state 
is flavor-polarized, a decrease in the activation would be 
expected — at least initially, before the condition for having a 
flavor-polarized ground state is violated. Here, the enhanced 
^-factor reduces the energy cost for slowly varying spin/flavor 
fluctuations according to Eq. ( 145 > . leading to an increase in 
optimal skyrmion size and hence a reduction in the optimal 
skyrmion energy. 



B. n = 1 Landau level 

Next, we will discuss quantum Hall ferromagnetism aiv = 
3,4, and 5, corresponding to a quarter-filled, half-filled, and 
three-quarter-filled n — 1 Landau level, (n — —1 is related by 
particle-hole symmetry.) We will again ignore Landau level 
mixing and project onto the n = 1 level. Since the wave 
functions now reside on both sublattices, we must return to 
the full interacting theory defined by Eqs. (Ill> and ( I12> . 

In the SU(4)-invariant limit with Hi = 0, the system will 
again be an exchange ferromagnet at filling factors v = i,A, 
and 5. It is interesting to explore the lowest-energy charge ex- 
citations here, and in particular to ask whether the transport 
gap is set by topological skyrmions or separated particle-hole 
excitations. In the case of Schrodinger equation Landau lev- 
els, it is known that skyrmions constitute the minimum-energy 
charge excitation only in the lowest Landau level. ^'' Since we 
are dealing here with two-component wave functions whose 
elements consist of lowest Landau level and next-lowest Lan- 
dau level Schrodinger equation wave functions, the answer in 
the present case is not a priori obvious. 

To address this issue, we calculate the energies of the gen- 
eralized SU(4) "spin-wave" branches in the SU(4)-invariant 
limit. While the convenient first-quantized formalism de- 
scribed in the previous subsection is no longer at our disposal, 
these energies can be computed directly using second quanti- 
zation. For simplicity, let us consider v ^ "i and assume the 
following symmetry-broken ground state. 



io)=n4 



R,m\ 



(49) 



where here c 



T-R,m 



adds an electron into the n = 1 Landau 
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(47) 
(48) 



level. (Our attention will be restricted to this filling factor 
since the energies are unchanged aiv = A and 5.) We evaluate 
explicitly the energy of spin wave excitations created by the 



operator S'-(q) = /^ e-^i-^[5°^(x) - iS'°^(x)], where 5°^ 
and S'^y are defined in Eq. j26> . This choice is arbitrary, since 
which branch we evaluate is immaterial due to the assumed 
SU(4) symmetry. Upon projecting onto the n = \ Landau 

level, we have S'^{q) -^ S'ct (q), with 



^^(q) 



-iq-x 



^[a>r^(x)]tci>;jr'W 



^ ['^J,K,;'^Tfl,m + C|i_;C|L,m]- 



(50) 



Here $p^j^ (x) are two-component wave functions defined in 
Eq. (I18> . with /, m labeling the angular momentum (we have 
suppressed the spin and flavor indices since the wave func- 
tions do not depend on these). A spin-wave in the projected 

space is thus given by |q) ~ Sa (q)|0), and has an excitation 
energy 



i?'(q) 



(q|7^o|q) (0|Ho|0) 

(o|o) ■ 



(q|q) 



(51) 



where TLq is the SU(4)-invariant part of the Hamiltonian pro- 
jected onto the n = 1 Landau level. Letting L? be the system 
size, we obtain 



Nik) = ^ 



V{k)sin^{i\khq/2)N{k), 



(52) 



X K='(y)]^'i>rHy) 



= (,-WBkf 



[l-{lBkf/^f 



(53) 



Note that replacing the wave functions in Eq. ( I53> by one- 
component Schrodinger equation Landau level wave func- 
tions reproduces the known spin-wave energies for the latter 
case.^^ Thus, Eqs. (|53 and Eq. ( I53> constitute a straightfor- 
ward generalization of the spin-wave energies to the case of 
two-component wave functions. 

Equation ( I52> yields the following energy dispersion at 
small momentum. 



7 [^ 



64 V 2 



Sc- 



(54) 
(55) 



This implies that a skyrmion-antiskyrmion pair costs an en- 
ergy 



^sk/ask — '^'^ ■ 



(56) 



The energy of a particle-hole excitation meanwhile is given 
by 

^;/. = ^oM = ^y|^c. (57) 

Comparing these, we see that E'^k/askl ^'p/h ^ '^/22 < 1; 
consequently, in SU(4)-invariant limit skyrmions will set the 
transport gap in the n — \ Landau level as well. 



1. Filling factor v = 3 

We will now incorporate the anisotropy terms in Hi to de- 
duce the probable ground states selected by interactions, con- 
sidering i/ = 3 first. Due to the Zeeman coupling, we will as- 
sume a spin-polarized ground state here. The ordering in fla- 
vor space remains to be determined. The most general flavor- 
polarized state favored by the long-range Coulomb repulsion 
can be expressed 



\9,^) = n 



771 — 



cos ^c\jj„^ + sm -e "Pcj^^ 



(58) 



The angles 9 and (p specify the polarization direction in flavor 
space. With the ground state spin-polarized, the only terms 
in Til remaining that break the flavor degeneracy are the vi 
and U2 interactions. Naively, the former favors a state with 
(Pstag) 7^ 0, while the latter favors having ( J±) ^ 0. It is easy 
to show, however, that both expectation values vanish for any 
choice of 6' and cj) in Eq. ( I58t . The vanishing of (pstag) follows 
from the fact that the n = 1 wave functions carry equal weight 
onbothsublattices. The vanishing of (J±) is a consequence of 
the fact that the elements (/)"^^ and x"^^ in Eq. ( II8I 1 are given 
by wave functions for different Schrodinger equation Landau 
levels; hence X]m['?^"r^]*Xrrr^ = 0- ^" other words, as op- 
posed to what we saw at i/ = — 1, the direct contributions from 
the interactions leave the flavor symmetry unbroken here. 

There is, nevertheless, subtler flavor symmetry -breaking by 
interactions that arises due to exchange. Specifically, vi fa- 
vors in-plane flavor polarization (6 — 7r/2), while U2 favors 
out-of-plane flavor polarization (e.g., 9 — 0). Denoting the 
projected vi and U2 interactions as 7Y„j and TYuj, respectively, 
we obtain the following variational energies: 



(7r/2,0|H.i+H„j7r/2,, 



PoL' 



U2P0 (59) 



{o,mv^+nu,\o,4') = ^uipo. (60) 

Since ui w \/3w2/4, it follows that the in-plane polarized 
state has lower energy so that the v = i ground state is given 
by 



1^=3)^ n [4/^,™+^*^^,™]!^"^)' (61) 

with a spontaneously chosen in-plane polarization angle (f). 
This state is characterized by a nonzero order parameter 
(■(/'^(r^ + iT^)-0) 7^ 0, and in general breaks translation, ro- 
tation, and reflection lattice symmetries for graphene. Using 
Eqs. (|5|i and (|6|l, one can show that this discrete symmetry 
breaking is due to lattice-scale fermion currents circulating 
around nearest-neighbor honeycomb plaquettes, with a spe- 
cific pattern determined by the polarization angle 0. 

Since this state spontaneously breaks U(l) flavor symmetry, 
a finite-temperature Kosterlitz-Thouless transition is expected 
here. The physics associated with this transition has been ex- 
amined in detail previously in the context of quantum Hall 
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bilayers, where the layer degree of freedom plays the role of 
flavor. -■'—'" One obstacle for observing such a transition in 
the bilayer problem is the finite tunneling amplitude between 
layers, which explicitly breaks the U(l) psuedospin symme- 
try that is analogous to the U(l) flavor symmetry in our case. 
Here, however, the U(l) flavor symmetry appears to be more 
robust. Graphene may therefore eventually provide a clean 
setting for studying the physics of a Kosterlitz-Thouless tran- 
sition in the quantum Hall effect. Further theoretical explo- 
ration of this phenomenon in the present context would cer- 
tainly be interesting. 



2. Filling factor u = 4 

Similarly to what we found at j^ = 0, the ground state 
at ly = 4 depends on the strength of Zeeman coupling and 
the on-site U relative to that of the interactions vi and U2- 
The former terms again favor a spin-polarized flavor-singlet 
state, while the latter interactions now favor an in-plane flavor- 
polarized spin-singlet. Evaluating the variational energies of 
these states, we find that the spin-polarized state has a lower 
energy if giisB + moPo > "iPo + M2Po/4. This inequality 
is satisfied even with a vanishing on-site U and a g-factor as 
small as unity. Thus we will assume here that v = 4 ground 
state is spin-polarized: 



=4>=n4H. 



m^'\L,m\ 



(62) 



We note, however, that the gaps for spin-waves out of this 
state are positive only \f gjisB + 2uopo > 2uipo + U2po/2, 
which provides a sUghtly more restrictive condition for having 
a fully spin-polarized ground state. With g — 1.4, the latter 
inequality is satisfied provided U > 2eV. 



3. Filling factor u = 5 

Finally, the situation at i^ = 5 is analogous to that atiy — S. 
The ground state will have all spin up states in the n — 1 Lan- 
dau level occupied to satisfy the Zeeman coupling. Due to the 
sublattice repulsion vi, the remaining spin down electrons in 
the ground state are expected to be in-plane flavor-polarized: 



5)^n[' 



'iR,m 



p«<*^t 



"i-L,m. 



1^-4). 



(63) 



This state exhibits the same type of lattice-scale order as the 
1^ — 3 ground state, and also spontaneously breaks U(l) flavor 
symmetry so that a Kosterlitz-Thouless transition is expected 
here as well. 



C. Experimental relevance? 

Initial sightings of the graphene quantum Hall effect ob- 
served integer quantum Hall states only at filling factors i^ = 
±2, ±6, ±10, etc.^ corresponding to filled nearly four-fold 



degenerate Landau levels. Quantum Hall ferromagnetism was 
thus entirely absent in these experiments. This is almost 
certainly due to a collapse of the exchange gap responsible 
for quantum Hall ferromagnetism by strong disorder, a phe- 
nomenon that has been established both theoretically '^'"'^iS 
and experimentally^i^. 

As discussed in the introduction, more recent experiments 
utilizing higher magnetic fields observed additional quantized 
Hall plateaus at i^ = 0, ±1, and ±4.^ A natural question to ask 
is whether the appearance of these plateaus is a manifestation 
of quantum Hall ferromagnetism. Nomura and MacDonald, 
who very recently derived an experimental criterion for real- 
izing quantum Hall ferromagnetism in graphene, suggest that 
this is indeed the case. However, the measured activation en- 
ergy ai V — ±4 was found to be dominated by the single- 
particle Zeeman splitting, suggesting that here the system is 
perhaps more appropriately characterized as aparamagnet. At 
present it is unclear whether the same is true for the plateaus 
at z/ = or 1/ = ±1. Measurements of the activation energies 
there will likely provide some hints as to their origin. 

In the following section, we will explore the possibility that 
the system is indeed paramagnetic at these newly observed 
quantized Hall plateaus, and that their origin is due to explicit 
symmetry-breaking terms in the Hamiltonian (rather than ex- 
change) as appears to be the case atiy — ±4. In an attempt to 
access this regime, we will assume that the effects of exchange 
are effectively "canceled" by disorder, and utilize Hartree the- 
ory to examine the spin- and flavor-splitting of the n = and 
n — 1 Landau levels by such terms. This analysis will hope- 
fully provide useful input for determining the true nature of 
the I' — Q and i' = ±1 states. 



V. PARAMAGNETIC REGIME 

It will be instructive to begin our analysis of the paramag- 
netic regime by applying Hartree theory initially in a clean 
system. Disorder effects will be discussed at the end of this 
section. Since in our treatment of the ferromagnetic regime we 
found that the U2 interaction dropped out at the n = Landau 
level and was effective only through exchange at n = 1, we 
will ignore this term here. Doing so will greatly streamline 
our discussion, and our results do not depend on this simplifi- 
cation. We retain the remaining symmetry-breaking terms in 
Til, and consider the following interacting Hamiltonian, 



n = d^x{ - ihy^^^lri'^D^ + if Dy\^j ~ gpeB ■ St 



+ Wo[pt,totPi,tot + PT:StagPi,i 



stag J 



uip: 



stag 



}> 



(64) 



where pa,tot and Pa,sta.g, are the total and staggered densities 
for spin a (e.g., Pt.stag = iplT^ri'^'ip^). We have ignored the 
finite range of the sublattice repulsion t^i, as a local form will 
now be adequate. Furthermore, we have dispensed with the 
manifestly SU(2)spin invariant form of the on-site U term in 
favor of a form which will be more convenient to work with 
here. Note also that we have dropped pieces of the on-site U 
interaction involving J±, since again these will be unimpor- 
tant. 
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Our aim will be to employ Hartree theory to obtain effective 
single-particle energy levels for the interacting system, and in 
particular to extract the resulting particle-hole excitation en- 
ergies which will be of potential relevance for experimentally 
measured activation energies. To proceed we decouple the 
interactions by linearizing in fluctuations of the total and stag- 
gered densities about their mean values, yielding an effective 
mean-field Hamiltonian 

-Hmf = I (f-^{-ihvij^[ri''D:,+f]yDy\iP 

~ gfJ-B^ ■ Stot + 2_^ [■^a,tPa,tot + ^a,sPa,sta.g\ 
— MsPstag) + ^^const- (65) 

The constant-energy term that arises upon decoupling is given 
by 



Er. 



= / d^x{ -■uo[(p|,tot)(Pi,tot) 

+ (PT,stag)(pi,stag)] +-Ui(pstag)^}. (66) 



Furthermore, the new coupling constants appearing in Eq. 
(165 > are defined as follows. 



A 



T/i:t 

Ms 



Uo{Pi/-\,tot) 

"o(Pi/T,stag) 

2ui(pstag)- 



(67) 
(68) 
(69) 



The Act^t and Xa,s terms represent effective chemical poten- 
tials for the total and staggered densities for spin a, while 
/is behaves as an overall staggered chemical potential. As- 
suming uniform expectation values for the total and staggered 
densities, it is now straightforward to diagonalize the mean- 
field Hamiltonian TYmf to obtain the desired effective single- 
particle levels for the interacting theory. Below we exam- 
ine the resulting level structure at integer filling factors in the 
n = Q and n = 1 Landau levels. 



A. n = Landau level 

The two-component n = Landau level wave functions 
for flavor A and spin a retain the same form as in the non- 
interacting case (see Sec. Illlt . Their energies, however, are 
now shifted as follows: 

^Vi.R = T 2-9^5^ + ^T/i,t--^T/i,s + Ms (70) 

^T/i,i = T^gP'BB + \yi^t + >^yi,s- Ps- (71) 

In particular, we see that the flavor degeneracy exhibited in the 
non-interacting case is generically lifted by interactions. Be- 
low we specialize to filling factors v = ±1 and i^ = to ob- 
tain the effective level structure in the respective ground states. 
We will find here that Hartree theory predicts the same ground 
states as in the ferromagnetic case, due to the fact that direct 
contributions from interactions rather than exchange selected 



the ordering in the latter case. Excitation energies, however, 
which are our main interest, will be dramatically modified by 
the loss of exchange. 



1. Filling factors v — ±1 

We will again consider only i/ = —1 since v = +1 is 
related by particle-hole symmetry. The mean-field energy at 
J/ = — 1 is minimized by filling all spin up, flavor R states (or 
all spin up, flavor L states). Setting (pt,tot) = -(pt^stag) = 
po and (px,tot> = (Pi,stag> == 0, we obtain the following ef- 
fective single-particle levels: 



E]R 
E\L 

EiR 



gpsB - 2uipo (72) 

gpsB + 2uipo (73) 

gpsB + 2uqpo - 2uipo (74) 

\gpBB + 2uipo. (75) 



By examining Eqs. ( I72> through (I75> . we see that the mini- 
mum energy required to make a particle-hole excitation out of 
the ground state is either 4wipo or giisB + 2uoPo, whichever 
is smaller. Again we note that from Eqs. Wil and M6\ we 
have uopo w 0.08C/[eV]B[T]K and uipo w 0.4B[T]K, 
where [/[eV] is expected to be on the order of a few electron- 
volts. Figure|2ja) illustrates schematically the effective single- 
particle levels in these two cases. 



2. Filling factor u — 

Just as we saw in the ferromagnetic regime, the v ^ Q 
ground state here depends on the strength of Zeeman cou- 
pling and the on-site U relative to the sublattice repulsion 
Ml. Evaluating the mean-field energetics for these states, 
we again find that the ground state will be spin-polarized if 
gpBB+2uopo > 4uipo, while a flavor-polarized ground state 
exhibiting lattice-scale CDW order occurs if gfiBB+2uQpo < 
AuiPq. 

Let us first deal with the case where the ground state is 
spin-polarized. Here, we set (p|,tot) — 2po and (px.tot) = 
(Pt/i-stag) = 0> yielding the following energy levels. 



psp 
^IR/L 



;gPBB 

-gfiBB + 2uopo. 



(76) 
(77) 



Since all spin up states are occupied, the energy for a particle- 
hole excitation is clearly given by gpsB + 2woPo- Con- 
sider alternatively the case where the ground state is flavor- 
polarized, and we have (pt/x,tot) = -(PT/i,stag) = Po cor- 
responding to occupying all flavor R states. The effective 
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FIG. 2: Schematic ejfective single-particle energy levels obtained 
within a Hartree analysis of the interacting theory at (a) // = —1, 
(b) u = 0, and (c) v = A. Here As = g^-sB + 2uQpo, while 
Auj = 4miPo- For v = —1 and 0, the levels on the left-side corre- 
spond to As > Auj , while the right-side corresponds to As < A„j . 



single-particles levels are now 



E. 



fp _ 



fp 

T/ii 



^^-gjisB + 2uqPo - 4uipo 



2' 
f 1 



(78) 
(79) 



Consequently, the minimum energy for a particle-hole excita- 
tion is SuiPq — 2uopo — gpsB. The single-particle levels for 
both scenarios are illustrated schematically in Fig.|3b). 



B. n = 1 Landau level 

We will make an additional simplification in our treatment 
of the n = \ Landau level and set A^g = /is = at the 
outset. This is not necessary and does not affect our results, 
and we do so only for the sake of simplicity. On physical 
grounds one expects these terms to vanish given that the non- 
interacting n = 1 wave functions carry equal weight on both 
sublattices. Furthermore, it is easy to show that in any case 
they do not lead to a splitting of the flavor degeneracy in the 
n = 1 Landau level. 

With this simplification, the wave functions are unchanged 
from the non-interacting case, while the energies for flavor A 
become 



Ef / I A 



1 



n/M -^2^^^BB + Xyl^t 



\/2ehv^B. 



(80) 



In contrast to the n = Landau level, the flavor degener- 
acy remains unbroken by interactions here. The important 



physical implication is that there are gapless charge excita- 
tions at 1/ — 3 and 5, and thus no quantum Hall effect would 
be predicted by Hartree theory at these filling factors. This is 
not too surprising since in our analysis of the ferromagnetic 
regime we found that interactions broke the flavor degeneracy 
at i^ = 3 and 5 only via exchange. 

The levels are spin-split, however, and there is a gap ?Av — 
4. All spin up orbitals are occupied in the ground state, so 
here we have (pi.tot) — 2po and (p|,tot) = 0. The effective 
single-particle energies are thus given by 

E\^ = -^g^BB + V2env^B (81) 

E[a - +-gtiBB + 2uoPo + V2dw^, (82) 

and the particle-hole excitation energy is gpsB + 2uoPo- The 
energy levels here are illustrated in Fig.|2jc). 



C. Discussion 

To summarize, our disorder-free Hartree analysis suggests 
that Zeeman coupling together with symmetry -breaking inter- 
actions can in principle induce quantum Hall states in a para- 
magnetic system at filling factors t/ = 0, ±1, and ±4, but not 
aiv = ±3 or ±5. The activation energies alv ^ ±4 obtained 
in Hartree theory — corresponding to one-half of the particle- 
hole energies — are 



A^=±4 = -gpsB + uqpq. 



(83) 



At filling factors v — ±.1 and v — Q, the activation ener- 
gies depend on whether (i) gpsB + 2uoPo > 4miPo or (ii) 
gpsB + 2uopo < 4uipo- In case (i) the activation energies 
at these filling factors are predicted to be 






2uipo, 
1 



V=o - 2^^^^ 



uopo 



while for case (ii) we have 



1 



-gpsB + uopo. 






(84) 
(85) 

(86) 
(87) 




We note that in the second case the activation energy at i/ 
is larger than that for ly ^ ±1 and ±4. 

These results will of course be modified by the inclusion 
of disorder. In particular, one may ask whether the lattice- 
scale character of the wave functions, which was essential for 
the flavor symmetry-breaking at n = and lack thereof at 
71 = 1, survives upon disordering the system. We expect that 
microscopic features of the wave functions will indeed remain 
robust, provided the length scale associated with the disorder 
is long compared with the lattice spacing. Furthermore, it is 
clear that assuming uniform expectation values for the total 
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and staggered densities as we did above will not be valid in 
the presence of disorder. This introduces a considerable de- 
gree of complexity, as the mean-field Hamiltonian in this case 
is not easily diagonalized and must be solved self-consistently. 
However, we expect that disorder can be effectively accounted 
for by a (possibly magnetic field-dependent) broadening of 
the effective single-particle levels found above, leading to a 
reduction in the activation energies we predicted by consider- 
ing a clean system. A more serious treatment of disorder and 
interactions will clearly be desirable, but will not be pursued 
here. 

Even at the crude level of our treatment, it is interesting to 
compare our predictions from Hartree theory with the high- 
field experiments mentioned above that resolved additional 
quantum Hall states at i^ = 0, ±1, and ±4. First, Hartree the- 
ory is consistent with the appearance of these quantum Hall 
states, and also can potentially provide a simple explanation 
for the conspicuous absence of those ?A v ~ ±3 and ±5; 
namely, a lack of flavor symmetry-breaking due to the form 
of the higher-Landau level wave functions. Second, the acti- 
vation energy predicted at z/ = ±4 is given by the Zeeman 
energy plus a contribution from the on-site U that is linear in 
the perpendicular magnetic field, the latter being unchanged 
by the addition of an in-plane magnetic field component [see 
Eq. (I83H . The measured energy gap for the spin-split Landau 
levels was similarly found to be given by the Zeeman energy, 
supplemented by a phenomenological Landau level broaden- 
ing varying linearly with the perpendicular magnetic field.- If 
one interprets this Landau level broadening as a contribution 
to the energy splitting arising from the on-site U, one obtains 
U ^ 2eV, which is of the right order of magnitude. Finally, 
if gfisB + 2uopo < Auipo, then one obtains activation en- 
ergies in Hartree theory which are identical at i/ = ±1 and 
±4, while a larger energy is predicted at i^ = [see Eqs. ( I83> . 
(I86> . and (I87H . It is intriguing to note that experimentally the 
quantized Hall plateau at i^ = sets in at B ~ IIT, while the 
i^ = ±1 and ±4 plateaus are resolved at similar fields around 
B w 17T, which may be consistent with the relative magni- 
tudes of the above Hartree energies. 

We emphasize that the preceding discussion is by no means 
intended to be conclusive, but only to suggest one possible 
mechanism for the appearance of the additional quantum Hall 
plateaus at high magnetic fields. Again, assessing the origin of 
the quantum Hall states at i/ = and v = ±.1 requires further 
experiments, which we hope this work may stimulate. Mea- 
suring the activation energies at these filling factors as a func- 
tion of an in-plane magnetic field may provide some guidance 
as to whether the quantum Hall states at these filling factors 



are due to quantum Hall ferromagnetism, explicit symmetry 
breaking, or perhaps some other mechanism. For instance, if 
gUsB + 2uoPq > -iuipo, our Hartree estimates predict that 
the i^ = gap will increase with an in-plane magnetic field 
component, while the gaps at i^ = ±1 will remain unchanged. 
On the other hand, for gpsB + 2uopQ < 4?iipo the gaps at 
both 1^ = and v = ±.1 vary with an increasing in-plane field, 
the former decreasing while the latter increases. 
VI. CONCLUDING REMARKS 

As technological progress enables the fabrication of higher- 
quality graphene samples, quantum Hall ferromagnetism will 
likely provide an interesting avenue of exploration along the 
road to the fractional quantum Hall effect. One of the remark- 
able aspects of quantum Hall ferromagnetism in graphene is 
the connection between quantum Hall physics operating on 
long length scales of order the magnetic length and lattice- 
scale physics. The clearest manifestation of this interplay oc- 
curs at filling factors v = ±.\ and possibly also v = 0, where 
lattice-scale charge density wave order, albeit weak, coexists 
with the integer quantum Hall effect. Establishing the pres- 
ence of such lattice-scale structure experimentally would be 
quite interesting, although doing so will likely prove chal- 
lenging as the signatures of such order are expected to be 
small. Another remarkable aspect of quantum Hall ferromag- 
netism in graphene that would be worth pursuing is the poten- 
tial for observing a Kosterlitz-Thouless transition, which we 
argued should occur at filling factors v — ±3 and ±5. Fur- 
ther experimental and theoretical studies of skyrmion physics 
in graphene would also be interesting. 

In the more immediate future, further experiments to deter- 
mine the origin of the additional integer quantum Hall states 
appearing at high magnetic fields would be extremely useful. 
We have provided here Hartree estimates for the activation 
energies in these states which, at least in principle, should be 
straightforward to compare with by exploring the dependence 
of the transport gap on an in-plane magnetic field. Finally, 
studies that incorporate disorder into our analysis of the para- 
magnetic regime in a more systematic way would of course be 
welcome. 
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